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Abstract 



We model the term structure of the forward defauh intensity and the defauh density by using Levy 
random fields, which allow us to consider the credit derivatives with an after-default recovery payment. 
As applications, we study the pricing of a defaultable bond and represent the pricing kernel as the unique 
solution of a parabolic integro-differential equation. Finally, we illustrate by numerical examples the 
impact of the contagious jump risks on the defaultable bond price in our model. 
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1 Introduction 

The term structure modelling in the interest rate and in the credit risk modelling has been widely adopted 
and extended since the original paper of Heath- Jarrow-Morton |[T6l . Notably, there have appeared many 
important papers (e.g. H] 01 El IH [HI [IH) incorporating jump diffusions to describe the family of bond 
prices or the forward curves as a generalization of the classic HJM model. 

In the credit risk modelling, the conditional survival probability associated to the default time is an 
important quantity for measuring default risk and studying valuation of credit derivatives. Let r be a 
nonnegative random variable defined on a complete probability space (17, T ^ P) equipped with a filtra- 
tion F = {Tt)t^Q satisfying the usual conditions. The conditional survival probability (CSP) is defined 
as StiO) = P(t > 6\J^t), t,9 ^ 0. To describe the term structure of the CSP, we can use both the density 
and the intensity point of view. On the one hand, as in El Karoui et al |[TOl . we assume that there exists a 
family of Tt (8) B(M+) -measurable functions {uj,9) — > at{uj,9) such that the CSP has the following additive 
representation: 



The family of random variables at( ) is called the conditional density of the default time r given Tt- On 
the other hand, similarly to the definition of forwai^d rate, we can use the "intensity" point of view and the 
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following multiplicative representation: 




■e 



St{9) = exp 



Xt{v)dv 



(1.2) 



where the J^t'^'S(M+) -measurable function (w, 9) — )• Aj (w, 6) is called the forward intensity. It is equivalent 
to assume the existence of the density or the intensity for all positive t and 0. We have the relationship: 



In the interest rate models, the time 9 is always larger than t and the forward rate has no economic inter- 
pretation for 9 < t. However, it is noted in ifTOl that to study what happens after a default event, we need 
the whole term structure of the conditional survival probability, that is, for all positive t and 9. One typical 
example is a defaultable bond where the recovery payment is effectuated at a given maturity later than the 
economic default date. 

In this paper, we consider the whole term structure modelling of CSP and the applications to the credit 
derivative pricing. In the credit risk models, the default contagion phenomenon is often modelled by positive 
jumps in the intensity process. We take this point into modelling consideration and propose a forward inten- 
sity driven by Levy random fields. In the existing Levy term structure models, in Filipovic et al. HBO, the 
authors consider forward curve evolutions as solutions of the infinite dimensional Musiela parametrization 
first-order hyperbolic stochastic differential equations driven by n-independent Levy processes or driven 
by a Wiener process together with an independent Poisson measure. In f9\, Eberlein and Raible present a 
class of bond price models that can be driven by a wide range of Levy processes with finite exponential mo- 
ments. This model was further applied to describe the defaultable Levy term structure and explore ratings 
and restructuring of the defaultable market. The driving process of the Levy term structure model in ^ was 
further extended to non-homogeneous Levy processes in HI. 

Motivated by those existing Levy term structure models and the random field models which are widely 
used to model various stochastic dynamics (e.g. |[Il|5l|6ll2lll5l[l7l[T8l), we suppose that the Levy random 
field in our model is a combination of a kernel-correlated Gaussian field and an independent (central) Pois- 
son random measure. The jump component described as Poisson measure is similar to that used in |[T2l . 
but it is not necessary to assume the exponential integrability condition for the characteristic measure under 
our framework (see Section 2). The kernel-correlated Gaussian field is more flexible compared to the Gaus- 
sian components without kernel-correlation considered in ifTSlfTTlfTSl . In fact, we can choose appropriate 
correlated-kernels of the Gaussian field so that the models considered in ||9l HT] [121 can be covered (see 
Remark 123] ). Note that it is not genetically tractable for pricing of defaultable bonds under infinite dimen- 
sional framework as in |[Tll.[T2l . Although we do not intend to consider the forward intensity under infinite 
dimensional framework as in IHIO, it has a close relationship between the (infinite dimensional) Wiener 
process and the kernel-correlated Gaussian field. Indeed, the kernel-correlated Gaussian field can product 
a cylindrical Wiener process by establishing appropriate Hilbert spaces (see Proposition 2.5 in |j7J). We 
deduce the dynamics of the CSP and the associated density in this setting. In particular, we emphasize on 
a martingale condition, which can be viewed as an analogue of the non-arbitrage condition in the classical 
HJM model. 

For the pricing of credit derivatives, we follow the standard general framework in Bielecki and Rutkowski 
|[3l . The global mai^ket information contains both the default information and the "default-free" market infor- 
mation represented by the filtration F, which is obtained by an enlargement of filtration. We are particularly 
interested in an economic default case, that is, the default does not lead to the total bankruptcy of the un- 
derlying firm and a partial recovery value is repaid at the maturity date of the bond in case of default prior 
to the maturity. To evaluate this "after-default" payment, we use the density approach in ifTOl and obtain 
that the key quantities for the pricing of a defaultable bond ai^e two pricing kernels, one depending on the 



at{9) = St{9)Xt{9). 



(1.3) 
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interest rate and the default density, and the other depending on additionally the recovery rate. We assume 
that both the short interest rate and the default density are modelled by the Levy random field model and 
are correlated between them. For the recovery rate, we analyze firstly the simple case where the recovery 
rate is deterministic and then the random recovery case. We show that the pricing kernel is related to the 
solution of a second-order parabolic integro-differential equation and we prove, based on a result of GaiToni 
and Menaldi llT4l . the existence and the uniqueness of the solution to the equation. 

The rest of the paper is organized as follows. We present our model setting in Section |2] and give the 
martingale condition. We then analyze the dynamics of the CSP and the conditional density in Section |3] 
In Section m we discuss the pricing of credit derivatives and in particular the defaul table zero-coupon bond. 
The two sections [5] and [6] focus on the pricing kernels. Finally, we present some numerical illustrations in 
the last section 17] 



2 Forward intensity driven by Levy random field 

In this paper, we adopt a random field point of view to model the forward intensity A((^) where both t and 
9 are positive. We consider a Levy random field on R+ x which is a combination of a Gaussian random 
field and a compensated Poisson random measure independent to Y'-'. Here denotes the time 
space and M"' is considered as a parameter space. 

We assume that the covariance of the Gaussian random field Y*-^ is given by a kernel measure c on M*^ 
which has a continuous and symmetric density on \ {0} with respect to the Lebesgue measure and such 
that c({0}) > 0. Namely for ((^i, <^2) G C^(^+ x 1^'^)^ 

Jo JR2d 

where by abuse of language c(^i — C2)d^id^2 denotes the measure on M'^ x M.'^ the inverse image of the 
measure c by the mapping from x M'' to M*^ which sends (^1,^2) to ^1 — ^2- The Gaussian random 
field defines a worthy martingale measure (see ||20l p.289] and E p. 190]). Let F'^ = (Jf )t5.o be the 
filtration satisfying the usual conditions which is generated by 

a{Y'^{[0,u] X A), u i:t, A e Bb{M'^)), t ^ 0, 

where Bb{^'^) denotes the set of all bounded Borel subsets of M'^. Let "P*^ be the predictable cr-algebra on 
Q X R_|_ associated to F*^ and be the linear space of all (g) i3(M'^)-measurable functions h such that 



JM2d 



1/2 

< +00 



for any T > 0. The stochastic integral h-Y^ is well defined for any h G <I>c- When c is the Dirac distribution 
concentrated on the origin, the stochastic integral: 

S(to,...,td) = y^([0,to]x---x[0,td]), (to,...,td) effif ^ (2.1) 

defines a (d + l)-parameter Brownian sheet. If in particular- d = 0, it becomes a standard Brownian motion. 

Denote the intensity measure of the compensated Poisson field Y^ by z^(d^)dt, (t, ^) G M+ x W^, where 
1/ is a (j-finite measure on W^. Let F^ = {J^^)t^o be the filtration satisfying the usual conditions generated 
by 

a{Y^{[Q,u]x A), u^t, A^B^iM!^)), t^O 
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and be the predictable cr-algebra on x associated to F^. Denote by "if^ the linear space of all 
(8) ;B(]R'^) -measurable functions such that 



T 

2, 



\g{t,0\'i^mdt 

JR'' 



< +00 



\\9\\u,T :=IE 

for any T > 0. The stochastic integral g ■ is well defined for any g G 

Let F = (/t)t>o be the natural filtration generated by the Levy random field, namely F := F*^ VF^. We 
describe the forward intensity by using the Levy random filed as the following additive HJM type model: 

dXt{9) = ^H{0)dt + [ at{e,OY'^idt,dO+ [ 7t-{e,0y''idt,d0, (2.2) 

where 

(1) fi = {nt{^); {t,e) e M^) is P ® -measurable and ^[\^it{e)\]dt < oo, where V is the 
predictable cr-algebra on x IR_|_ associated to the filtration F, 

(2) a = {at{e,^); {t,9,() e M+ x M+ x R'^) is ® x M'^) -measurable and for any 6 € M+, 

(3) 7 = {jt{e, C)> 0; {t, e, C) G M+ X M+ X R'^) is «) B{R+ x M'^)-measurable and for any 9 € M+, 



The model (12.21 ) can also be written in the integral form as 



Xt{e) = Xo{9)+ l' fis{9)ds + 





f [ as{e,Oy''{ds,dO+ f f ls-{e,Oy''ids,dO (2.3) 
Jo Jo Jw 

where both stochastic integrals with respect to and are F-martingales with mean zero, Ao(-) is a 
deterministic Borel function on IR+. 

Similarly to the classical HJM model, in the above Levy field model (12.21 ). there exists a relationship 
between the drift coefficient /x and the diffusion coefficients a and 7 due to the fact that, for any 9^0, the 
conditional survival probability process {St{9) = exp ( — Xt{v)dv) , t ^ O) should be an F-martingale. 
We call this relationship the martingale condition (MC). Let us introduce the following notation: 

I^{t,9):= ^it{v)dv, I^it,9,0-= Mv,Odv, and I^it,9,0-= lt{v,C)dv, 
Jo Jo Jo 

where (t, 9, ^) e M+ x M+ x R'^. 

Theorem 2.1 For all 9 ^ and T > 0, one has 

[ E[\I^{t,9)\]dt<oo, I^i-,9,-) e<^c, Lf{;9r) e^,, and e'^-^-'^'-) - 1 G (2.4) 
Jo 

Moreover, the process family [St{9) = exp ( — Xt{v)dv^ , t ^ O) is a family of ¥ -martingales if and 
only if the following condition is satisfied: 

(MC) V O 0, tit{9) = f c7t{9, Ci)Ut, 9, 6)0(6 - 6)d6d6 

(2.5) 

7i(0,e)(l-e-^-(*''^'«V(d6- 
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Proof. The proofs for the first three statements in (12.41) are similar. We only provide the details for the third 
one. For any T > 0, we have 



T r rT 

E|/^(t,e,^)|2i/(dOdt = 
Jr'^ Jo 

T 



E 



E 



JR'i 



'yt{vi,S,)dvi / 'yt{v2,(,)dv2 



1 

< - 
- 2 



^0 Jo 

T 



^0 



which is finite since 7. (u, •) G for any v ^ 0. For the last assertion in (12.41) . note that 7t(^, and 
thus 



< \Ij{t,0,O\ 



for all {t, e, G X R+ X M'^. 

We now prove that the condition (MC) is equivalent to the martingale condition for {St{0), t ^ 0). In 

fact 



St-{e) 

1 



{t,d)dt- [ i^{t,e,OY^{dt,dO 

JRd 



+ 7T / iAt,e,^i)Ut,e,C2)c{Ci-C2)d^id^2dt 
+ / (e-^^(*"'^'«) - i)y^(dt,dO 
+ f (e-^.(*.^'?)-l + /^(t,0,e))z.(dOdt, 

JR'i 

so the martingale condition of {St{9), t ^ 0) is thus equivalent to the following equality 

i^{t, 9) = \ [ h{t, 9, 9, 6)c(a - e2)d6d6 

2 jR2d 

which is equivalent to (MC). 



(2.6) 



□ 



Remark 2.2 Consider the particular case where c? = 0, c is the Dirac measure, and u = The condition 
(MC) becomes 



V^^O, ^it{9) = <Jt{e) / <Jt{v)dv. 

Jo 

This corresponds to the non-arbitrage condition in the classical HJM model where the forward intensity is 
driven by a standard Brownian motion. 
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Remark 2.3 There exist random field models in the literature. We make below some comparisons. The 
forward intensity model (2.2) can be extended to the following form: 

dXt{e) = fiti9)dt+ [ at{e,C)Y'^{dt,dC)+ [ jt^{e,OY''idt,dO 

+ / ^tA0,OiY''idt,dC) + uidOdt), 

where a. (9, •) € <I>c, 7.(6', Oljl Ki} ^^d 7.(6', Oljl-^i} ^ ^u, for each fixed 9^0. Under the model dTT] ). 
the con^esponding martingale condition (MC) will be changed accordingly. We next consider a special 
form of the predictable random filed with separable variables: 

where (a* (6*); {t, 9) G M^) is a real-valued predictable random field, (74(6*) = {^j (9),... , ^ti^)); (t, &) G 
M^) is a R'l^-valued predictable field and (f){() is a deterministic measurable function on W^. In this case, 
the extended model ( 12.71 ) can be rewritten as 

dXt{9) = {fit{0)-a)dt + at{9)Y^{dt,^{*)) + {^t{9),dLt), (2.8) 

where a G M, (•, •) denotes the inner-product on and 

dLt = adt+ [ CY^idt, + / ^{Y^idt, + i^idOdt) 

is a non-Gaussian Levy process if the characteristic measure is a Levy measure. If (p = 1, then y"^(]l[o x 

(/){*)) becomes a Brownian motion when the correlated-kemel is Dirac. Choose appropriate smooth function 
(j) as in the proof of Proposition 2.5 in [71, then y'^(l[o t] x <^(*)) becomes a cylindrical Wiener process. 
Thus we recover the Levy interest rate term structure models considered in ||9j[TTl[l2l, if the Levy measure 
u satisfies the exponential integrability condition. We next give a comparison of our Levy random field 
+ Y^ introduced previously in this section with existing Levy fields in literature. 

1. As in ( 12.11 ). the field Y'^ + Y^ can be reduced to a Brownian sheet in Walsh ||20l . when the kernel c 
is Dirac and the characteristic measure u = (hence Y^ = 0); 

2. the field Y^ + Y^ becomes a so-called "colored" space-time white noise model established by IH, 
when the kernel c(^) = |^|~" with < a < d and v = 0; 

3. the fractional space-time white noise (fractional in space and time in white) used in |fT9T| corresponds 
to the field Y'^ + Y^ with the kernel c(C) = h{2h - with i < /i< 1, d = 1 and z> = 0; 

4. the Poisson sheet in m corresponds to the field Y^ + Y^ + v{d^)dt with c = and z/(^) = z6i [d^] 
where z > is single point and 5i is the Dirac measure concentrated at 1. The Gamma sheet in lUl 
is the field Y^ + YP + iy{dC)dt with c = and u{dS,) = ^l^^^oyzd^ where d = 1 and z > is a 
single point. 

3 Conditional survival probability and density 

In this section, we concentrate on the conditional survival probability {St{9), t ^ 0) and the conditional 
density {at{9), t ^ 0). Here we specify a cadlag version of the martingale {St{9), t ^ 0) for any 9^0. 
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In fact, to show that the integral Xt{v)dv defines a cadlag process, we need a stronger assumption on the 
process X{6) in order to apply Lebesgue's theorem. The cadlag version of S{6), if well defined, should have 
a universal version of its predictable projection as follows: 



Xt-iv)dv]. 



Thus 



St{9) := lira exp (- [ X„.{v)dv] (3.1) 
qm+ \ Jo J 



git 

defines a universal cadlag version of the martingale S{9). 

We observe from the equality (12.61) that, under the condition (MC), the conditional survival probability 
admits the following dynamics: 



dStje) 



-[ Ut,e,C)Y''{dt,dO+ [ (e-^-(*"'^'«)-l)y^(dt,de), (3.2) 

where Sq{9) = exp(— XQ{v)dv). 

For ^ 0, we denote by M{6) the martingale defined as 

dMtie) = - [ h{t, e, i)Y^{dt, di) + [ (e-^-(*-'^'«) - i)y^(dt, dc), Mo{e) = o. (3.3) 

With this notation, S{6) / So{6) is the Doleans-Dade exponential of the martingale M{6). Moreover, denote 
by m{9) the martingale defined by the dynamics: 

dmti9) = - [ ati9,0Y''idt,dC)- [ 7*- Oe-'^^*~'''^^>^''(dt, d^, moW = 0. (3.4) 
Observe that the following relation holds 

re 

Mt{e) = I mt{u)du. 



Jo 

We then consider the dynamics of the conditional density of default given in (11.11 ). 
Proposition 3.1 The conditional density process a{9) admits the dynamics: 

datiO) = at^{e)dMt{e) - St^{e)dmt{e) (3.5) 

or equivalently 

oit-{0) Xt-{0) 



Proof. Keep the martingale condition (MC) in mind. The dynamics (13.51 ) is derived by employing Ito's 
formula to a{9) = X{6)S{9) for each positive 9 fixed. □ 

An important property in the credit analysis is the immersion property, or the so called (H)-hypothesis, 
which means that an F-martingale remains a G-martingale. The (H)-hypothesis is satisfied if and only if 
at{9) = aQ{9) or equivalently Xt{9) = X0{9) for any t ^ 9. In the random field setting, by (12. 31 ). this is 
equivalently to 

/ / as{9,OY''ids,dO= [ I 7s-(^,0^''(ds,de) = 



7 



for t ^ 9,OT equivalently 

at{e,O = and ^t{e,O = z.(de)-a.e. 

for t > 6. Note that the martingale condition (MC) then implies that = for t > 0. 

We recall that the F-intensity process A of the default time r coincides with the diagonal forward inten- 
sity, i.e. At = Xt{t). It is closely related to the Azema supermartingale: 

St = St{t)=F{T>t\J^t), 
which is also called the survival process of r. 

Proposition 3.2 Let M be the W-martingale having the dynamics 

dMt = - [ Ut,t,OY^{dt,dO+ [ (e-^-(*-'*'«)-l)y^(dt,de). 

Jwi jRd. V / 

Then ^ 
where £{M) is the Doleans-Dade exponential of M. 

Proof. The Azema supermartingale 5 has a multiplicative decomposition of the form St = Lt exp(— Xsds) 
(see ifTOl Proposition 4.1]), where L is an F-martingale having the following dynamics 

dLt = exp ( [ Xsds] dLt, 

with 

ft 







Lt = - at{u) - au{u)du. 
Jo 

By Proposition 13.11 together with ( 13.31 ) and (13.41 ). 

dLt = - [ f A{t,e,OdeY''{dt,d^)- [ f B{t,9,^)ddY''{dt,d^), 
Jm.''- Jo Jm.''- Jo 

where 

A{t, 6, = -at^ {o)h{t, e, c) + St^{e)at{e, 0, 

By integration by part, we obtain 

- [' A{t,9,C)d9 = -St-{t)Ut,t,C), 
Jo 

- ['B{t,e,Ode = 5t_(t)(e-^- - 1). 
Jo 

Moreover, the Doob-Meyer decomposition of S is given by 

St = l + Lt - / au{u)du, 
Jo 

which implies that 

dLt dLt dSt 

1— = = + -^id*- 

Jt- ^t- 

By (13.11) . one has 5j_ = St^{t). Hence the martingale L is the Doleans-Dade exponential of M and the 
assertion follows. □ 
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4 The pricing of defaultable bonds 



In this section, we focus on the pricing of credit derivatives. In general, a credit sensitive contingent claim 
can be represented by a triplet (C, G, R) (see Bielecki and Rutkowski IH) where the J^r-measurable ran- 
dom variable Ct represents the maturity payment if no default occurs before the maturity T, and G is an 
F-adapted continuous process of finite variation such that Go = and represents the coupon payment. 
Differently from the case where the default payment occurs at r immediately, we assume that in the eco- 
nomic default case, the default (or the recovery) payment takes place, after a period of legal proceedings, 
at the maturity date T later than the economic default date r and admits the form Rt{t) where Rt{-) is 
Ft ® S(]R+)-measurable. 

The global market information is described by the filtration G = {Qt)t^Q, Gt = TV cr(r A t), which is 
made to satisfy the usual conditions. The value at time t ^ T of the contingent claim (G, G, R) is given by 
the following -conditional expectation: 



Vt = E( 



{CTl{r>T} + liryuye-^^^'^^'dG, + l|.^r}^r (r)) e" 



(4.1) 



where Q denotes a risk-neutral pricing probability measure and the interest rate r = (r^; t ^ 0) is an F- 
adapted process. The following result computes Vt using -conditional expectations. The first two terms 
result from 131 and the third one from ifTOl . With an abuse of notation, we denote in the following the F- 
conditional density of r under the risk-neutral probability Q by {at{0),t ^ 0). The general result on the 
density under a change of probability measure is given in ifOl Theorem 6.1]. 



Proposition 4.1 We suppose that the economic default time r admits a conditional density w.r.t. the filtra- 
tion F, denoted by at{-) under the risk-neutral probability measure Q. Then the value of the credit sensitive 
contingent claim (G, G, R) is given by 



Vt 



{CtSt+ I RT{u)aT{u)du)B^^ + [ SuB-^dGu 
Jt Jt 



Tt 



+ l{r^t}BtE^ 



^TiO)-^jJ-Bj,^\Tt 
at{6) 



(4.2) 



where St = Q(r > t\Tt) = at{9)d6 and Bt = exp(/Q r^ds). 

Proof. The -measurable random variable Vt can be decomposed in two parts 14 = ^{T>t}V t+^{T<t}Vt{T) 
where Vt is -measurable and Vt{-) is T (8 S(M+) -measurable. On the set {r > t}, we use Jeulin-Yor's 
lemma (see Q) and the conditional density to obtain 



Vt 



1 

S't 

1 

s't 



{CtSt+ r RT{0)aT{e)de)e- + S^e' ''^'^'dG^ 
Jt Jt 



Tt 



On the set {r ^ t}, by HOl Thm 3. 1], we have 

Vtir) = E 
which complete the proof. 



at(d) 



T 



r.ds 



Tt 



□ 
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Note that for the pricing of the two "before default" payment terms (C, G), the quantity 



— = exp 



and hence the intensity A play an important role. However, for the default recovery payment R (which 
depends on r), the "after-default" density at{9) where t ^ ^ is needed. This point has been discussed in 
lITOl . In the following of this paper, we adopt the density approach for both the bef ore-default and after- 
default pricing. 

We consider in particular a defaultable zero-coupon bond of maturity T with C = 1 and G = 0. Its 
price at t ^ T is given by 



'Mt>t} + '^{t^t}Rt{t)) exp ^- ^ r,ds^ 



(4.3) 



The following result is a direct consequence of the previous proposition. We first introduce the following 
price kernels: 



K2{t,9) 



at{e) 



aT{0) exp 
RT{O)aT{0) exp 



T 



rgds 



(4.4) 
(4.5) 



where t ^ T, and 6^ ^ 0. 



Corollary 4.2 Using the conditional density of t under Q, the price (14.31 ) of the defaultable zero-coupon 
bond at time t has the following representation: 



P{t,T) 



1 



{r>t} 



T 



Ki{t,e)de + j K2{t, 



MO) 

St 



de 



l{r^t}K2it,T). (4.6) 



We will identify the above price kernels in the next two sections with different settings. 



5 The first pricing kernel 

In this section, we study in detail the pricing kernels ( 14.41 ) and ( 14.51 ) when the random interest rate is de- 
scribed as an extended Vasicek model. We suppose in this section the after-default recovery payment is 
deterministic. The case where the after-default recovery payment is random will be considered in the next 
section. 

Firstly we recall the forward intensity model (12.21) and assume that the F-predictable random fields 
(/I, <T, 7) are deterministic in (12.2b . We then express the instantaneous interest rate process r = (r^, t > 0) 
as the following extended Vasicek model under the risk-neutral pricing measure Q : 

drt = K{5-rt)dt+ f pt{Oy^{dt,dO+ [ MOY''{dt,dO, (5.1) 

where k > 0, 6 > 0, and p.(-) and (/).{■) are deterministic volatiUty functions, assumed to belong to $c and 
"ifi^ respectively. In the particular case where d = 0, 4>t{C) = the volatility function p.(-) = p > is 
constant, the interest rate r satisfies the classical Brownian-driven Vasicek model: 

drt = k{6 - rt)dt + pdWt, (5.2) 
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where is a standard Brownian motion. 

Similarly to the solution form of the Ornstein-Uhlenbeck stochastic differential equation, the extended 
Vasicek model (15.11 ) also admits an explicit expression as follows: 



n = roe 



Jo Jr<i 



+ 



(5.3) 



10 Jm.<i 

where ro > denotes the deterministic initial interest rate value. 

Next we compute the first pricing kernel in (14.41) . For ^ 0, we introduce the following integro- 
differential operator Aq acting on functions with three variables t, x and y which ai^e differential in t and 
second-order differentiable in {x, y) : 

dK dK d'^K 

AeK{t, X, y) = n{5t{9) - x) — {t, x, y) + a{t, 9) — {t, x, y) + aii(t)— ^(i, x, y) 

ox oy ax^ 



(5.4) 



where 



a{t,9)= fiti9)- 



+ / [K{t,x + M0,y + 7t{9,0)- K{t,x,y) 

dK dK 1 

-</'t(e)7^(i,a;,y) -7t(^,e)i^(t,2;,2/) i^idO, 
ox oy J 



ati9,c)i4t,e,0ciC-0dCdC 



7i-(0,O(l-e-'-(*''''«V(de), 



a22it,0) 



ai(0,ei)cTt(e,6)c(6-6)deid6, 



ai2{t, 9)= [ at{9, 6)^(6)0(6 - 6)deid6. 



Remark 5.1 Recall the martingale condition (MC) given by (12.51) which has been assumed throughout the 
paper. We have the coefficient a{t, 0) = for the partial derivative ^ under (MC). 

We introduce the following assumptions where we have fixed 6^0. 
Assumption 5.2 (1) There exists q E (0, 1) such that 

(i) the functions aii(-), a22{-, 9) and ai2(-, 9) are |-Lipschitz on [0, T], 

(ii) there exists a Borel function Jq on (which could depend on 9) such that 

max{|</.i(0 - 0.(01, \lt{9,i)-^s{9,0\} < MOlt - sl"^' 

and 



/ 



1 + MO 



i^(dO < 00. 
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(2) and |7t(^, ^)| ai^e uniformly bounded from above by a Borel function Jo(0 such that 

(3) There exists a constant f3{9) > such that, for any (x, y) E and any t £ [0, T], one has 

au{t)x'^ + 2ai2(t, e)xy + a22{t, 9)y^ ^ /3{e){x^ + y^). 

Then we have the main result of this section. 

Theorem 5.3 Let 9^0 be fixed. Under Assumption 15.21 the Cauchy problem 
dK 

— {t,x,y)-xK{t,x,y) + AeK{t,x,y) = 0, K{T,x,y)=y (5.5) 

has a unique solution K, where the integro-differential operator Aq is defined in ( 15.41 ). Moreover, the 
following equality holds 



aTi9) exp ^- Vgds^ 



T 

' 'Tt =St{9)K{t,rt,Xt{9)), (5.6) 



where St{9) = Q(r > is CSP and Xti9) is the corresponding forward intensity under the pricing 

measure Q. 

Proof. The first assertion comes from a general result of Garroni and Menaldi lfT4l Theorem II.3.1]. Let o 
be as in Assumption 15.21 We shall actually prove that the Cauchy problem (15.51 ) with a terminal conditiorlj 
K{T, ■,■) = ip e C''(M?) has a unique solution in the Holder space C^+2'2+'?([o, T] x M^) by constructing 
a contractible operator. The case of (15.5b with unbounded terminal function ip{t, x,y) = y will be treated 
by taking hmits. We recaU that C^+2'2+9([o, T] x M^) denotes the vector subspace oE C^'2([0, T] x M^) 
of functions / such that 

:= ||/||i,2 + E idid^M)t,^,i,+a+b+2c-i) + E < +00, 

l<a+b+2c<2 a+b+2c=2 

where for any function g : [0, T] x — M and any f3 € (0,1), 

/\ \9{s,z) - g(t,z)\ sr^ \g{t,z) - g(t,w)\ 

{g)t,f} ■■= sup sup -g , {g){x,y),(3 ■■= sup . 



'The expression ^(R^) denotes the vector space of all bounded functions / on which are Holder continuous of order q 
(namely, such that ||/||sup + ||/||g < +00), where 



II, := sup 



\h{z)-h{w)\ 



z,»eH- \z-w\i 



"The expression (7^'^([0, T] x M^) denotes the vector space of all continuous functions / on [0, T] x such that 

II/II1.2 := Yl Wdtd^d'yflU.^ < +00. 

a+b+2c<2 

The vector space C^''^{[0, T] x R^) together with the norm || • ||i,2 form a Banach space. 
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The vector space C^+3'2+9Qo, T] x M^) together with the norm 
Let I0 be the integro-differential operator defined as 



ll+f ,2+f/ 



form a Banach space. 



{IeK){t,x,y) 



K{t, x + MO,y + ItiO, 0) - K{t, X, y) - MO^ - ItiO, 0^ 

ox oy 



For K £ j^] ^ M^) and £ C"?(M2), let e^{K) be the unique solution of the Cauchy problem 



dF 

— -xF + Ae{F) = Ie{K), F{T,x,y) = i;{x,y), 



(5.7) 



where Ag denotes the differential operator 

92 



dx"^ 



aii(i)^ + «i2(i, + a22{t, 0)—^ + K{5t{9) - x)— + a{t, 6) — . 



d_ 

' dx 



d_ 

' dy' 



Denote by C2''i'([0, T] x E^) the vector space of functions / on [0, T] x such that 



/||sup+ sup \\f{-,x,y)\\i+ sup •, 

{x,y)&? te[0,T] 



< +00 



which is a Banach space with respect to the norm 



Since K G C^+3'2+'?([o, T] 



), by the 



AssumptionE2(l.ii)and(2), we obtain that /^(K) G C2'9([0, T] xM^) (see niXemmaII.1.5]). Therefore 



the existence and uniqueness of the solution 0^(K) G C^"'"2'2+'?([0, T] x M^) to ( 15.41 ) comes from the 
classical theory of parabolic partial differential equations (e.g. US). Moreover, the solution verifies the 
following Holder estimate ( |[T4l Theorem 1.2.1]) 

IIG^(i^i) - e^(i^2)||i+a,2+, Ci\\Ig{K, - 7^2)11,,, (5.8) 
which holds for all Ki,K2£ C^+2'2+'?([0, T] x M^) such that Ki{T, ■, •) = K2(T, •, •) = ip, where Ci is a 

l^) (see HH Lemma 



constant independent of ip. 

For arbitrary e > 0, the following estimate holds for any K G C^"''2'2+'?([o, T] x 
II. 1.5]) 



|/eWIL „ ^ e||V^,,j^)i^|U^^ +C(e)(||ir||,,, + |lV(,,,)(i^)|| 



(5.9) 



2,9 " ^-^''y/ "J, 9 

where the constant C(e) only depends on e. Denote by the subset of functions in C^^2'2+'3Qq^ j"] x M^) 
whose restriction on {T} x coincides with V'- By choosing e > small enough, we obtain from (15.81) 
and ( 15.91 ) that is a contracting operator on the complete metric space C^,, provided that T is sufficiently 
small. Hence for sufficiently small T, the operator 0^ has a unique fixed point and therefore the Cauchy 
problem 

dK 

-^(i, X, y) - xK{t, X, y) + AeK{t, x, y) = 0, K{T, x, y) = y) 

has a unique solution. For general T, it suffices to divide [0, T] into a finite union of small intervals and 
resolve the Cauchy problem progressively. 

For the terminal function ip{x, y) = y, we can take, for each integer n > 1, a function ^„ G C^(]R2) 
which coincides with (/? on the ball B^. of radius n centered at (0, 0). For any n ^ 1, let be the unique 
solution of the equation ( 15.51 ) with terminal condition Kn{T,x,y) = ipn{x,y). By a maximum principle 
for the equation (15.51) (see Ill4i Theorem II.2.15]), for n ^ m, coincides with Km on the ball Bm.. By 
taking K = Kn on [0, T] x Bn we obtain a global solution to the Cauchy problem (15.5b . The uniqueness of 
K also results from the maximum principle. 
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We now prove the second assertion that K satisfies the equality (I5.6I ). To this end, we compute the 
denominator of the pricing kernel (14.41) by introducing a change of probability measure: 



By B ayes' formula and (11.3b . we have 

fT 



axiO) exp 



_ St{0) 

Soiey 



(5.10) 



Xrid) exp 



rAs 



Note that, by Girsanov's theorem (see H Theorem 3.3]), under the probabihty measure Q^, 



defines a Gaussian field with correlated kernel c on R'^, and 

y^(dt,dO := l^^(dt,de) + (1 - e-^^(*'^'^V(de)dt 

defines a compensated Poisson random measure with predictable compensator e~^T^^*'^'^^i^(d^)dt. Then the 
dynamics (15.11 ) of the interest rate r can be rewritten as 

drt = Ki$t{e) - n)At + / ptim'^iA^. de) + / m)y''{'^^. do, 

where 

\(e) = b^K-^ [ pt(e)/.(i,0,e)c(c-e)dcde + A^-i / 0t(e)(e-^-(*'''«)-i)Kdo, 

and the dynamics (12.21 ) of the forward intensity rate can be rewritten as 
where 

{e) = f,t{0)- [ at{e,c)Ut,e,Oc{c-mdC- [ ^t-{e,m-e-'-^''''^^HdO- 



Note that the forward intensity process A(^) is a (Q^, F)-maitingale for each 6 fixed. Assume that 



Xt{0) exp 



r,ds 



K{t,rt,Xt{e)), 



where the function K(t, x, y) is sufficiently regular. Then Ito's formula applied to the (Q^, F)-martingale 



exp 



rAs]K{t,rt,Xt{e)) 
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yields 



-rtK{t,rt,Xtm + ^ii,rtAtie)) + K{U0)-rt)^it,rt,Xtm+M0)^{t,rt^ 
+ au{t)^{t, rt, Xt{9)) + a22(t, 0)^{t, n, Xt{e)) + a^it, 0)^^{t, n, Xt{e)) 



+ 



- MO^it, rt, XtiO)) - jt{0, O^it, rt, Xt{e))] u{dO = 0. 



Conversely, if K is the solution to 



dK 

— -:,K + AeK = 0, K{T,x,y) 



then one has 

fT 



Er 



ariO) exp ^- r^ds^ 



Sti9)K{t,rt,Xtie)). 



Thus we complete the proof of the theorem. □ 



Accordingly the pricing kernels (14.41 ) and (14.51 ) are given by 

Ki{t,e) = ^K{t,rt,Xt{e)) (5.11) 

K2{t,e) = ^RT{9)K{t,rt,Xt{e)) (5.12) 

where t ^ T and ^ ^ 0. By Corollary I4.2[ we obtain immediately the pricing formula for the defaultable 
zero-coupon bond. 

Remark 5.4 Concerning the pricing kernel at the left side of the equality (15.61 ). one possible alternative 
way is to solve it directly by using the dynamics of the density at{6). However, in view of (13.5b . the 
con^esponding solution K{t,rt, St{6), Xt{9)) will include three variables apart from time variable. The 
main advantage of the change of probability method (15.101 ) is that we obtain the solution function in the 
form K{t,rt, St{9), Xt{6)) = St{6)K{t,rt, Xt{6)). This indeed decreases the dimension of variables for 
our pricing kernel function and is important in the numerical computation. 

Remark 5.5 If the interest rate r is independent of the forward intensity, hence independent of the density, 
then the computation of the pricing kernels is easier. Denote by B(t,T) the price of the standard zero- 
coupon bond, i.e. B{t,T) = EQ[exp(- // r^ds)|7't]. Recall that we have assumed the recovery rate 
deterministic in this section. Then 

Ki(t,0) = ^ii^l|M), K2{t,e) = RTi9)Bit,T) 



which imphes that the time-t value (14.61 ) of defaultable zero-coupon bond has the following representation: 

This quantity serves to measure the default risk including both the default probability and the loss given 
default. We also notice in (15.131 ) that the recovery corresponds to a "recovery of face value" since it can be 
written as the quotient between the defaultable bond and an equivalent default-free bond. 
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Remark 5,6 The zero coupon price B{t,T) := 



exp 



K{t, rt) where K{-,-) is the unique solution to the following integro-differential equation: 



can be given in the form 



dK 



dK 



d^K 



dt 



dx 



dx"^ 



dK 



xK+ — + K{5-x)— + aii{t)—^+ / K{t,x + 4>t{0)-K{t,x)-(t)t{i)^{t,x) u{di) = Q 



dx 



with the terminal condition K{T, x) = 1. If there is no jumps in rt, i.e., (ptH) = 0, then the above equation 
becomes 

dK ,dK ^^d'^K ^ 

-xK + — + k{S- x)— + aii(t)^^ = 0. 



Its unique solution is 

K{t, x) = exp 



dt 

1 _ e-«(^-*) 



dx"^ 



{5-x)-5{T-t)+ / au{u) 



du 



where aii(t) is given in (15.41 ). Thus we obtain the following equality B{t, T) = K(t, rt), which is similar 
to the classical case. 



6 Random recovery rate and the second pricing kernel 



In this section, we consider the general case for the pricing kernel (14.51 ). where the after-default recovery 
payment is random as an extension to the previous section. 

Bakshi et al. lH assumed that the recovery rate is related to the underlying intensity as the following 
form: Rt = wq + wie~^^, wq,wi ^ {), wq + wi ^ 1 and A is the intensity process of default. In a similar 
manner, we assume that Rt{0) is of the form 

RT{e) =wo + u;ie-^(^^(^)) , ^ ^ (6. 1) 

where \t{G) is the forward intensity implied by (11.21) under the pricing measure Q, wq, wi satisfy the same 
condition as above and / is a non-negative function which is locally Holder continuous of positive order. 



Proposition 6,1 Let 9^0 be fixed. Under the Assumption \5.2\ the pricing kernel (14.51 ) is given by 



K2{t,e) 



-k{t, rt, XtiO)) + ^Kit, rt, Xti9)), 



Xt{9) ^ ' " - - ' Xt{e)' 
where K and K are respectively solutions to the partial integro-differential equation: 



(6.2) 



dK 
~dt 



{t, X, y) - xK{t, X, y) + AgK{t, x,y) = 



(6.3) 



under the terminal conditions K{T, x,y) = y and K{T, x, y) = ye 

Proof. Similarly to Theorem 15.51 the equation ( 16.31 ) with the terminal condition K{T, x, y) = ye 
admits a unique solution K. Moreover, by a change of probability measure we obtain 



'fiy) 



oiT{9)e 



-/(At(9)) 



exp 



Tt 



St{e)K{t,rt,\t{e)). 



Hence the formula (16.21 ) follows from the following relation (see (14.41) . (14.51) and (16.11 )) : 



K2{t,e) = ^,K,{t,e) + ^^ 



at{e) ' ^ ' at{9) 
where Ki{t, 6) is the first price kernel (14.4b . 



□ 
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Corollary 6.2 Under the Assumption \5 .2\ the price of the defaultable zero-coupon bond is given by 



P{t,T) =l|,>t} ' *^ ^-'^ - ^ ' ' ' 



k{t,rt,Xt{0))d9+ / -^(woK{t,rt,Xt{e)) + wiK{t,rt,Xt{e))y9 



T i^t Jt 'Jt 

+ ^\r^t} -T-i-T woK{t, n, Xtir)) + wiK{t, n, At (r)) 
where K and K are given in Proposition 16. W respectively. 

7 Numerical illustrations 

In this section, we illustrate our previous results by numerical examples. We are particularly interested in the 
contagion phenomenon. More precisely, we shall analyze in detail the the jump part in the default density 
dynamics and its impact on the defaultable bond pricing. 

In the numerical example, we consider the dynamics of the default density described by (13.51 ) and we let 
the martingale m{6)he given by 

dmt{e) = -(Jt{e)dWt + / 7t-(^)ee-«ioS*-Wd-y^'(dt,dO, mo(0) = 0, (7.1) 



with W = {Wt, i ^ 0) being a standard Brownian motion independent of the Poisson measure . 
Compared with (13.41 ). the corrected kernel c of the Gaussian field is the Dirac measure and d = 1, the 
volatihty coefficient (Tj (0, ^) = (Tt{6) does not depend on ^ and the jump amplitude coefficient is given by 

7t(0, ^) = 7t (^)Cl{5>o} where 74(6*) > 0. Recall in addition that Mt{9) = Jq mt{u)du and 

dat{9) = at^{9)dMt{9) - St^{9)dmt{9). 

To illustrate the impact of the jump part on the defaultable bond price P{t, T) given by (14. 61 ). we first 
consider the case when the martingale m{9) has no jumps, i.e., 7 = 0. We then include the jump part in 
the density dynamics. We use the initial default density given by qo(^) = Ae"'^^ with A being a positive 
constant. 

In the coming tests, we suppose that cJt{9) and ^t{9) are deterministic and we use the following forms 
of the coefficients and the characteristic measure in (17.1b . 



' at{9) = a{9-t)+, a > 
< ^t{9) = b{9 - t)+, b > 0, 
, = 4e-«/-%>o}de, C>0,zu>0. 

We assume that both the recovery rate i? G [0,1] and the interest rate r are constants and define B{t,T) = 
e~''(^~*) for ^ t ^ T. By Remark [531 the defaultable bond price P{t, T) given by ( 14.61 ) admits an explicit 
form. Since the quotient P{t,T)/B{t,T) equals the constant R on the set {r ^ t} in this case, we only 
study the pre-default part on {r > t} in (15.131) . which is denoted by P{t, T) henceforth and is given by 

P{t,T) = S(t,r) (l-(l-i?)i|^^^^y (7.2) 
V h M9)d9 J 

The main task is then to approximate the integral at{9)d9 by a finite sum '^^J:f/^_^i ^ * c^tii * A). 
Here we choose A = 1/100 and = 10/A. We perform 10^ experiments to compute the J'f -measurable 
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random variable P{t,T). In each experiment, we first generate the underlying Brownian motion and the 
central compound Poisson process. Then for each 6 G {zA; i = 1,2, ... , N/ A}, we compute at^{6) on 
{ti = iAf, i = l,2,..., t/At} with At = 1/100. 
The preferred parameter values are as follows: 

t = 0.5, T=l, r = 0.05, R = 0.4, b = l, ( = 10, X = 0.1. 

Figure [Uplots the kernel estimations of the densities of P{t, T) given by 

1 ^ 

fpix):=-Y,fh{x-m,T)), (7.3) 

i=l 

where Pi{t, T) is the price obtained in the i-th experiment, fh{x) = exp f— , and h = 1.06sfc/c^^/^ 
is the bandwidth, with being the sample standard deviation. From Figure [T] we find that the existence 
of the jump risk will increases the decentrality of the price. The right tail of the price distribution becomes 
fatter and fatter as the mean jump size w increasing. 



Figure 1: 
0.001. 



The (normal) kernel estimations of the price densities for tu = 0, 0.0002, 0.0006, 0.001, 0.002 and a = 
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Figure 2: Pe(0.5, 1) := E[P(0.5, 1)] as a function of w with A = 0.01, 0.03, 0.1, 0.3 and ct = 0.001. 
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Figure |2] shows the mean of the price Pe(0.5, 1) := E[P(0.5, 1)] as a function of zu for different values 
of A. We observe that the defaultable bond price is a decreasing function of the intensity A, and also of the 
mean jump size w. Hence, when there is larger default risk of the underlying asset itself (with larger A), the 
corresponding bond price is smaller. Furthermore, when there is more significant counterparty risks, that 
is, when there is a larger contagious jump in the density (larger w), then the bond price will also decrease. 
Both observations correspond to the reality on the market. 

Figure [3] shows the mean of the price Pe{t, 1) := K[P{t, 1)] as a function of t. It is noted that the 
numerical illustration of the quantity P{t,T)/B{t,T) discussed in Remark 1531 is very similar to that of 
P{t, T), since B{t, T) here is a deterministic function B{t^ T) = 6"*^^^^*^ which is close to 1. We observe 
similar results as in the previous test: the counterparty jump risks in the density will decrease the bond 
prices. 

Figure3: Pe(t,l) :=E[P(i,l)] as a function of t Right hand side is the relative price Pe(t, 1) := E[P(f, 1). 




0.2 0.4 0.6 0.8 1 ■ 0.2 0.4 0.6 0.8 1 

t t 

In the last graph, we show the quoted bond price at the initial time t = as a function of the maturity 
time T for different values of intensities. Again we observe that the bond price is decreasing when there is 
larger default risks and for long term bonds. 

Figure 4: P(0, T) as a function of T. 
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